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WEIGHT STUDY OF PARTIALLY SEGMENTED DIRECT-CONDENSING RADIATORS 

FOR LARGE SPACE POWER SYSTEMS 
by Roy L. Johnsen 
Lewis Research Center 

SUMMARY 

An analytical investigation has been performed to provide an insight into the weight 
increase of a direct -condensing radiator for a potassium Rankine cycle space power 
system as the system size was increased from 1 to 10 megawatts electric while holding 
constant the component efficiencies and all cycle temperatures and pressures. 

The largest single contributor to the radiator weight is the armor necessary to pro- 
tect fluid passages from meteoroid puncture. The two most important factors affecting 
armor thickness are the meteoroid population estimate and the radiator design survival 
probability. Design survival probabilities of 0. 9, 0. 95, and 0. 99 were used with a high 
meteoroid population estimate, Whipple’s (1961), and a low estimate, Watson's (1956), 
for 1, 10, 20, and 100 radiator segments. 

The radiator weight is sharply influenced by the meteoroid flux estimate and the de- 
sign survival probability. For the unsegmented case, using Watson’s (1956) flux esti- 
mate, the radiator weight ranged from 1. 5 to 3 pounds per kilowatt electric for the 1- to 
10-megawatt power range for a 0. 9 radiator design survival probability, and 2. 5 to 
6 pounds per kilowatt for a 0. 99 design survival probability. When Whipple's (1961) flux 
estimate was used, the radiator weight ranged from 7 to 17 pounds per kilowatt for a 
0. 9 design survival probability and from 16 to 42 pounds per kilowatt for a 0. 99 design 
survival probability . 

Segmentation of the radiator provided significant weight savings over the unseg- 
mented radiator, especially when high survival probabilities were demanded. Generally 
the weight advantages of segmentation increased with increases in system power level, 
design survival probability, number of segments, and severity of the meteoroid flux es- 
timate. When Whipple's (1961) meteoroid flux estimate and a 0. 99 design survival prob- 
ability for a 10-megawatt system were used, the radiator with 100 segments was about 
one-fourth the weight of an unsegmented radiator. 

With segmentation the ratio N /N of the design value of surviving number of seg- 

s 



ments to the total number of segments that yielded minimum weight was found to be 0. 8 
or higher. In this study, however, the manifolds were not segmented. Had the manifold 
been segmented, N /N would tend to be lower and conceivably might approach the 0. 75 
previously reported in the literature. 


INTRODUCTION 

The heat rejection aspects of electric generating systems for space application have 
received considerable attention in recent years. It is recognized that a large fraction of 
the power system weight would be taken up by the waste heat radiator if the conventional 
fin and tube type were used. Several radiator studies have been reported in the literature 
for systems up to several megawatts of generated power. A partial list includes refer- 
ences 1 to 6. It would be difficult to relate the results of these studies to one another 
since each study used different temperature levels, materials, and criteria for calcula- 
ting protective armor thicknesses. 

The analytical study presented herein is an attempt to determine weight trends that 
might be expected for both segmented and unsegmented direct-condensing fin and tube 
radiators for a Rankine cycle turboelectric system for the range of power levels from 1 
to 10 megawatts electric. No allowance was made for the weight of any necessary radia- 
tor support structure. The working fluid was potassium and the materials used were 
beryllium for the fins and tube armor and columbium for the tube liner. When high and 
low meteoroid flux estimates were used, weights obtained when segmentation of the tube 
panels was employed were compared to the unsegmented radiator weights. Of interest 
also, was the optimum ratio of design value of surviving segments to the total number of 
segments. Reference 7 had determined the optimum ratio of surviving to total number of 
segments for an armored radiator without fins and manifolds. 

The study objective was to glean trends from many detailed radiator designs that 
could be used to obtain rough estimates of the effect of power level and segmentation on 
radiator weight without having to make detailed calculations. Therefore, for compara- 
tive purposes, identical cycle conditions were maintained throughout the power range 
considered. 


PROCEDURES AND ASSUMPTIONS 

The procedure and the equations used are covered in some detail in appendix B. The 
development of the equations used to calculate such items as fin thickness, fin length, 
manifold weight and surface area, tube and manifold pressure drops, and fin and tube 
angle factors are shown in appendixes C to H. 
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It was assumed that saturated vapor at 1960° F entered the turbine, the turbine ef- 
ficiency was 0.8, the generator efficiency was 0. 9, and a turbine exit and condensing 
temperature of 1400° F had the corresponding vapor pressure of about 15 psia. The tur- 
bine temperature ratio was taken from reference 8 for minimum radiator area for the 
1400° F condensing temperature. Selection of cycle temperatures and efficiencies of the 
turbine and generator dictates the turbine exit quality, fluid mass flow, and radiator heat 
load. The calculated turbine exit quality was 0.85, the fluid mass flow was 6.88 pounds 
per second per megawatt electric, and the radiator heat load, including 100° F subcool- 

P7 

ing, was 1.8817x10 Btu’s per hour per megawatt electric. 

The following is a brief discussion of other assumptions and restrictions that in- 
fluence the study results. 


Materials 

The only materials considered were beryllium for the fins and meteoroid protective 
armor, and a 0. 03-inch columbium liner for the fluid passages. 


Radiator Panel Configuration 

Two configurations were used, both radiating from both sides. One was a two-panel 
configuration (fig. 1(a)) with a vapor manifold at one end and a liquid manifold at the 
other end. The other configuration (fig. 1(b)) with four panels had a central vapor mani- 
fold and liquid manifolds at the outer extremities. All manifolds were tapered. 





(a) Two-panel configuration. 


(b) Four-panel configuration. 


Figure 1. - Radiator panel arrangement. 
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Angle Factors 


In the radiant heat -transfer treatment, inter radiation and interreflection between 
surfaces were not considered. The angle factors for both the fin and tube took into ac- 
count only the direct radiation to space. The error introduced by such an assumption has 
been discussed in reference 9, and the maximum error for the results reported herein 
was approximately 8 percent. That is, the calculated radiant heat transferred from the 
fins and tubes was less than the actual amount by a maximum of 8 percent. 


Fins and Tubes 

Only constant thickness fins and constant diameter tubes were considered even 
though tapered fins and tubes may have offered some savings in weight. Reference 10 
shows a maximum weight savings of 10 percent with tapered fins, but did not include 
manifold weights. Had the manifold weights been included, the weight savings would 
have been even less. 


Factors Affecting Armor Thickness 

Radiator weight is closely tied to the amount of armor necessary to protect against 
meteoroid damage. The uncertainty of the numbers of particles that would be encountered 
and of the damage that would be inflicted by any one meteoroid particle makes a calcula- 
tion of armor thickness somewhat speculative. Several estimates of the meteoroid popu- 
lation are shown graphically in reference 11 and differ by orders of magnitude. The es- 
timates for the density of meteoroids vary from 0. 05 to 7.9 grams per cubic centimeter 
(ref. 12). Due to these wide differences in meteoroid flux and mass estimates, two sets 
of meteoroid data from reference 11 were used in order to encompass a range. Arbi- 
trarily selected were Whipple’s (1961) and Watson's (1956) meteoroid flux estimates. The 
meteoroid flux and mass estimates and the other important inputs used in the calculation 
of armor thickness are the following: 

(1) Meteoroid flux. 

(a) High estimate, Whipple's (1961): F-. = 9.27xl0“ 7 m" 10 / 9 particles/(ft 2 )(yr) 

ir 

(eq. for curve in ref. 11). The meteoroid density recommended for use with this 
flux was Pp = 2.7 grams per cubic centimeter. 

(b) Low estimate, Watson's (1956): = 1. 19x10"^ m”* particles/ (ft 2 ) (yr) (eq. 

for curve in ref. 11). A meteoroid density of p^ = 0. 5 gram per cubic centimeter 
was arbitrarily assumed. The armor thickness is only weakly influenced by particle 


4 



density. The meteoroid flux used is the dominant factor in influencing the calculated 
armor thickness deemed necessary. 

(2) Penetration criterion. The model proposed in reference 11 was used. The equa- 
tion is shown as equation (B7 ) of appendix B . 

(3) Survival probability, 0.9, 0.95, and 0.99. 

(4) Mission duration, 2 years. 


Manifold Fluid Velocity 

Tapered manifolds maintained constant fluid velocities. The fluid velocity was 
chosen, then the manifold internal diameter was calculated to maintain the velocity 
picked. For most of the study the fluid velocities used were 200 feet per second in the 
vapor manifold and 2 feet per second in the liquid. If the sum of the calculated manifold 
pressure drops exceeded 0.25 psi, the velocities were reduced to 150 feet per second 
for the vapor and 1 foot per second for the liquid. If the sum of the pressure drops still 
exceeded 0. 25 psi, the vapor velocity was further reduced to 100 feet per second. 


Other Assumptions 

Radiant input. - The effective sink temperature was assumed to be 0° R. 

Working fluid temperat ure. - The condensing working fluid was assumed to be iso- 
thermal everywhere. The fluid stagnation temperature was used to obtain the fluid ther- 
modynamic properties. 

Internal tube wall temperature. - The tube internal wall temperature was assumed 
to be equal to the fluid stagnation temperature. 

Subc ooling . - The condensed liquid was subcooled 100° F. 

Emissivity. - The external radiating surface emissivity was e = 0.9. 

PARAMETERS FIXED AS A RESULT OF PRELIMINARY CALCULATIONS 

Preliminary calculations were made in an attempt to reduce the number of variables 
to be used throughout the study. A discussion of these preliminary findings and reasons 
for picking the values selected follows. 


Generalized Length Parameter 


In reference 13 the value for the generalized fin length parameter 


L that yielded 
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Figure 2. - Effect of generalized length 
parameter on radiator weight and area. 
Whipple's flux (1961); system power out- 
put, 10 megawatts; probability of re- 
ceiving no meteoroid puncture during 
mission duration, 0.9; individual tube 
length, 40 feet; tube static pressure 
drop, 1 psi; two-panel configuration. 
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Figure 3. - Specific weight and ratio of fin span to tube outside diameter as 
function of number of tubes. Watson's flux (1956); system power output, 

1 megawatt electric; probability of receiving no meteoroid puncture during 
mission duration, 0.9; individual tube length, 15 feet; tube static pres- 
sure drop, 1 psi; generalized length parameter, 1; two-pane! configuration. 


the minimum volume fin was approximately 0. 92 for the case where the sink temperature 
was 0° R. The value of that would yield the lightest weight radiator when tubes and 
manifolds were included was not obvious. Values for L of 0.7, 0. 92, 1. 0, and 1. 1 

B 

were used to make some preliminary radiator weight and area calculations. For each 
value of Lg the number of tubes was varied until a minimized weight and corresponding 
tube diameter were obtained, while everything else such as individual tube length, radi- 
ator heat load, and tube pressure drop were held constant. Figure 2 is typical of the re- 
sult of plotting the percent deviation from minimum weight against L . For these pre- 

O 

liminary calculations, the minimum weight usually occurred at a value of L close to 

5 

unity. Figures 11 and 12 of reference 10 showed combined fin and tube heat rejection 
per unit weight, plotted against the ratio of half fin length to tube outer radius, to be 
maximum for a conductance parameter of about 1. The conductance parameter of refer- 
ence 10 is identical to L 2 /e. Using an emissivity e of 0. 9 the conductance parameter 

o / 
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of unity (of ref. 10) corresponds to a value for of about 0. 95. In reference 10, the 
particular value of ^fi n/ /( D 0 /2) at which the maximum heat rejection occurred was differ- 
ent for each case shown. 

While the approach taken in the study reported herein was somewhat different from 
that of reference 10, changing the number of fixed length tubes, while maintaining a con- 
stant Lg, changes the percentages of the total heat load rejected by the tubes and by the 
fins. This in essence is the same as changing the ratio 2jf^ n /D Q . This can be seen 
from figure 3 where specific weight and 2JZ^ n /D Q are plotted, for one example, against 
the number of fixed length tubes. Decreasing N^. increases and D Q but 2£^ n /D Q 
also increases. Thus for finding minimum radiator weight (constant L ), changing N. 
is equivalent to changing 2i^ n /D o . All of the subsequent radiator weights were obtained 
by using a value of L = 1. While a value close to unity provided the minimum weight, 
it should not be overlooked that by using a lower value of the projected radiator area 
may be decreased somewhat for only a slight radiator weight penalty. For example, re- 
turning to figure 2, when a value of = 0. 7 is used, the area could be reduced by per- 
haps 9 percent with a weight increase of about 2 percent. For some applications it might 
be desirable to accept the higher weight to get the smaller area. 


Pressure Drop 

The pressure drop across the radiator affects its weight in two ways. The size of 
the fluid passages is a function of the pressure drop, a high pressure drop permitting 
small passages. High pressure drop, however, reduces the average condensing temper- 
ature as well as adversely affecting the weight of the whole power system by depressing 
the cycle efficiency; this effect has been ignored in this analysis. Furthermore, a high 
radiator pressure drop results in a low condensate pump inlet pressure which may make 
necessary additional subcooling to prevent pump cavitation problems. Figure 4 shows a 
weight proportionality factor, derived by dividing radiator weight by the weight of an 
arbitrarily chosen radiator, for three tube pressure drops plotted against the number of 
tubes for one power level as an illustrative example. The minimum weights for each 
tube length were about the same for the 1- and 2-psi tube static pressure drops and 
higher for the 5-psi pressure drop. From figure 4 a tube static pressure drop of 1 psi 
was chosen for the remainder of the study. This may not be the optimum pressure drop, 
but the trends being sought should not be influenced significantly. Another 0. 25 psi max- 
imum was allowed for the sum of the pressure drops in the vapor and liquid manifold. 
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Figure 4. - Effect of pressure drop, tube length, and number of tubes on 
weight. Meteoroid flux, system power output, probability of receiving 
meteoroid puncture during mission duration, generalized length param- 
eter, and number of panels are constant. 


DISCUSSION OF RESULTS 

Radiator calculations were made at four power levels: 1, 2, 5, and 10 megawatts 
electric. For each power level several tube lengths were considered, and for each tube 
length specific radiator weight was plotted against the number of tubes . This was done 
for the unsegmented radiator for all combinations of the three design survival probabil- 
ities, two meteoroid flux estimates, and the two panel configurations considered. These 
plots are too numerous to include, but figure 5 is included as a representative example. 
Tube length and the number of tubes were varied until the combination that yielded mini- 
mum weight was obtained for each power level. In some instances the minimum weight 
for a particular tube length could not be obtained because the calculations were stopped 
when the fluid velocity at the tube inlet reached 1300 feet per second. Reference 14 in- 
dicates that the vapor velocity for critical two-phase flow (steam water mixtures) is 
lower than the sonic velocity for the vapor. The equilibrium sonic velocity for potassium 
saturated vapor (ref. 15) for the 1400° F condensing temperature is about 1450 feet per 
second. The 1300 feet per second limit was chosen arbitrarily to stay below the vapor 
equilibrium sonic velocity by a safe margin. 
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Figure 5. - Effect of tube length and number of tubes on radiator specific weight. 

Whipple’s flux (1961); unsegmented; probability of receiving no meteoroid punc- 
ture during mission duration, 0.95; tube static pressure drop, 1 psi ; gener- 
alized length parameter, 1.0; two-panel configuration. 

Effect of Number of Tubes and Tube Length on Tube Internal Diameter 

The internal diameter of the tube for any particular radiator may be of interest be- 
cause of its influence on the area vulnerable to meteoroid puncture as well as heat- 
transfer considerations. As previously mentioned, a constant tube static pressure drop 
of 1±0. 05 psi was used for the entire study. The internal tube diameter and approximate 
fluid velocity at tube inlet for any radiator shown in figure 5 can be obtained from fig- 
ure 6 (p. 10). It is seen from figure 5 that for any given power level increasing the tube 
length decreases the number of tubes at which the minimum weight occurs and that indi- 
vidual internal tube diameter increases (see fig. 6). Figure 6 also shows that for a 
given tube length the internal tube diameter is relatively insensitive to changes in the 
number of tubes, even though the fluid velocity at tube inlet varies considerably. This 
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Figure 6. - Internal tube diameter as function of tube length and number of tubes. Tube static pressure drop, 1 psi. 

indicates that as inlet velocity increases, increases in friction pressure drop are offset 
by increases in the momentum pressure rise. 


Planform Geometry 

The ratio of the width of the radiator to the length in the vehicle axial direction may 
be important from a vehicle integration standpoint. Flexibility in choosing either dimen- 
sion without having to pay a severe weight penalty would be desirable. Plots of specific 
weight against the ratio of total width to radiator length (fig. 7) indicate that some lati- 
tude in the choice of planform layout does exist. The minimum weight points from plots 
similar to figure 5 were plotted to show specific weight as a function of the ratio of radi- 
ator width to length in figure 7. Each curve represents the locus of the minimum weight 
points . 

For some of the curves the specific weight was not very sensitive to the width to 
length ratio. To cite one example, in figure 7(a) for the four -panel configuration for 
2 megawatts electric and 0. 9 design survival probability using Whipple's (1961) meteo- 
roid flux estimate, the width to length ratio ranged from 0.4 to 10 with a specific weight 
variation of only about 10 percent. Generally it can be said that as the survival probabil- 
ity is increased or the meteoroid flux estimate is increased, increasing the meteoroid 
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Figure 7. - Effect of ratio of total width to length on radiator specific weight. 
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protection problem, sensitivity to changes in width to length ratio decreases and the 
minimum weight panel tends to approach a square. 


Influence of Power on Weight 

The minimum weights for each combination of power level, survival probability, 
meteoroid flux, and number of panels were taken from plots similar to figure 5 (p. 9) 
and used as the basis of the discussion in the section that follows. 

When the minimum radiator weights at each of the four system power levels were 
plotted against the system electrical power, the points on a log-log plot fell in a straight 
line, which suggested an expression of the form 

WccP^ (1) 




System power output, P e , MWe 

(a) Two-panel configuration. (b) Four-panel configuration. 

Figure 8. - Effect of meteoroid flux estimate, survival probability, and system power level on radiator specific weight 
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as used in reference 5. (All symbols are defined in appendix A. ) If the value of the con- 
stant exponent b is known, radiator weights can be scaled from one power level to 
another with some degree of confidence. 

Equation (1) was rearranged to 



pb-1 

e 


( 2 ) 


which is the form of the curves plotted in figure 8 for the unsegmented radiator. The use 
of Whipple’s (1961) meteoroid flux estimate, for the 0. 9 radiator survival probability, 
yields radiator specific weights of 7 pounds per kilowatt electric for a 1 -megawatt system 
up to 18 pounds per kilowatt electric for a 10-megawatt system. For the 0. 9 survival 
probability, using Watson’s (1956) meteoroid flux estimate, the radiator specific weight 
ranges from about 1.3 pounds mass per kilowatt electric for a 1 -megawatt system to 
about 3 pounds mass per kilowatt electric for a 10-megawatt system. For the higher 
0. 99 survival probability, the specific weights were heavier than the specific weights as- 
sociated with the 0. 9 survival probability by a factor of two or more. 

The value of b in equation (2) was approximately 1.4 when Whipple's (1961) meteo- 
roid flux estimate was used, and only slightly less when Watson's (1956) meteoroid flux 
estimate was used. Increasing the design survival probability from 0. 9 to 0. 99 increased 
the value of the exponent b only about 2 percent. Also, the value of b was essentially 
the same for both the two- and four -panel configurations. The value of the exponent re- 
ported in reference 5 was 1 . 6 for a beryllium radiator for a heat load between 1 and 
10 megawatts thermal. Results plotted in reference 6 would lead to a value of less than 
1. 2 for a direct condensing radiator and slightly higher for a liquid filled radiator, both 
using copper as the fin material. 


Segmenting 

Reference 7 indicates that segmenting the radiator to reduce the amount of necessary 
protective armor should offer significant weight savings. Affecting the amount of pro- 
tective armor necessary for a given overall survival probability are the number of seg- 
ments N and the design value for the number of segments N remaining unpunctured at 
the end of the mission. Probabilities S of 0. 9, 0. 95, and 0. 99 that at least N g out of 
N segments would not be punctured for the 2-year mission duration were assumed to 
calculate the armor thickness for varying N while N was held constant at 10, 20 and 
finally 100 segments. It should be pointed out that the manifolds were not segmented, 
only the tube banks, and that calculations were made for the two-panel configuration only. 
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The equation for calculating the armor thickness is shown in appendix B. Furthermore, 
it was stipulated that all of the radiator heat load was to be rejected from N_ segments 
and the manifolds. 

Reference 7 showed a figure of merit, which is proportional to weight, plotted as a 
function of the ratio of N to N from which it was concluded that to obtain minimum 
radiator weight the optimum ratio of N g to N was approximately 3/4. Reference 7 did 
not include manifolds and fins. An attempt was made to compare some of the results of 
reference 7 with the results of this study which does include manifolds and fins. 

In reference 7 it was assumed that the required armor thickness t was proportional 
to the cube root of meteoroid particle mass m^, which in turn was inversely proportional 
to the meteoroid flux F^: 


t a 00 3 cc FC; 1//3 (3) 

The F % is the meteoroid flux having mass m or greater . The armor thickness t 
is calculated such that particles smaller than will not penetrate . Therefore, is 
also the average puncture rate per unit of vulnerable area. Use of the general case 
where 


F 


> 


°c m 


-0 

P 


(4) 


would have led to 


t oc m V 3 oc FC 1 / 3 ^ (5) 

a, p 

Thus reference 7 has used the special case where 0 = 1. For use where 0/1, the 
equations in reference 7 were rewritten to include 0. The figure of merit r? is then de- 
fined as 



( 6 ) 


The relation between weight and rj then becomes 


W 


N a 
s v 


— « 0(N s a v t) 1/3/3 


(7) 
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where W’ is the weight of the segmented prime surface radiator. Reference 7 suggests 
s 

that for analysis of a given vehicle and mission N g a v r may be considered to be con- 
stant, which makes 77 a very useful figure of merit in that it is directly proportional to 
weight. Carrying this further, the ratio of segmented to unsegmented weight can be ob- 
tained by 


W' 

s 77 


W u 


% 


( 8 ) 


where 7 ^ is the value obtained from equation (6) when the radiator is unsegmented, so 

that N and N are both unity. The values for 77 presented in reference 7 were for the 
s 

special case where j3 = 1 and this is the value of the exponent used in Watson’s (1956) 
meteoroid flux estimate. Whipple's (1961) meteoroid flux estimate has a value of 
£ = 10/9, therefore 77 was recalculated using the new value of /3; the W g /W^ was also 
calculated. 

The ratio of segmented to unsegmented weight W /W , which includes fins and un- 
segmented manifolds obtained by actually calculating the radiator weights for the seg- 
mented and unsegmented radiator is plotted in figure 9 (p. 16) for the high flux estimate 
and in figure 10 (p. 17) for the low flux estimate. Also plotted in figures 9 and 10 is the 
ratio W’/^, for the prime surface radiator, obtained from equation (8). Just as the 
survival probability P(0) of the unsegmented radiator was set at 0. 9, 0. 95, and then 
0. 99, the overall probability S that N g out of N segments as well as the manifolds 
would not be punctured was set at 0. 9, 0. 95, and then 0. 99. Expressed mathmatically for 
N identical and independent segments, 


N 


S = P(0) m x P(0) N = P(0) m x 


N 


' m 


y - 

/—! n! 

n=N 


NJ 

(N - n)! 


p n d - p> N - n 


(9) 


The probability that the manifold would not be punctured P(0) m and the probability that 
N of N segments would remain unpunctured P(0) M were equally weighted. With S 

S JN g 

assigned, P(0) and P(0) N are simply the square root of S. The individual segment 
. iN s 

probability of nonpuncture p, which when used in equation (9) yielded the desired value 

of S, was used when calculating tube armor thickness. 

Two generalizations can be made about all of the plots in figures 9 and 10. First, 

the ratio N /N at which the minimum weight occurs is somewhat higher than the ex- 
s 

pected 0.75 from reference 7. When using Whipple’s (1961) meteoroid flux estimate, the 
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Figure 9. - Effect of number of segments surviving and total number of segments on segmented radiator relative weight. Whipple’s flux 
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Figure 10. - Effect of number of segments surviving and total number of segments on segmented radiator relative weight. Watson's flux 




minimum weight occurred at N /N of 0.8 or higher, and when using Watson's (1956) 

s 

meteoroid flux estimate, the minimum occurred at N /N closer to 0. 9. This shift to 
higher values is probably due to the fact that the manifolds were not segmented. The 
initial assumption that the entire heat load be rejected from N g segments means that 
N - N g segments are redundant at the start of the mission. Increasing the ratio N g /N 
from 0. 75 to 0. 8 or 0. 9 reduces the number of redundant segments, which consequently 
reduces manifold length and weight. 

The second generalization that can be made is that for any given survival probability, 
number of segments, N /N ratio, and meteoroid flux estimate the value of W /W de- 
creases as the system power level increases; that is, segmenting is more advantageous 
at higher power levels. 

When figure 9 and Whipple's (1961) flux estimate are used, it can be seen that the 
plot of W* /W' is a good indicator of the relative weight for 10 and 20 segments, but for 
100 segments the level of W„/W is higher than W'/W'. The curve for W'/W' ap- 
plies to a radiator where 100 percent of the weight is protective armor. When Whipple's 
(1961) meteoroid flux estimate is used for the radiator with fins and unsegmented mani- 
folds, the tube armor comprises only about 60 percent of the total radiator weight. As 
the number of segments is increased, the tube armor comprises a smaller percentage of 
the total weight, which makes W^/W^ less applicable as a weight trend indicator . When 
Watson’s (1956) meteoroid flux estimate is used (which is considerably lower than 
Whipple's, 1961), the tube armor comprises a lesser percentage of the total weight of the 
radiator with fins and unsegmented manifolds. Thus W'/W' as a weight trend indicator 
is even less applicable when Watson's (1956) flux estimate is used for the radiator with 
unsegmented manifolds as can be seen in figure 10. 

It is seen by examining figure 9 that the relative weight W /W for 10, 20, and 
100 segments with Whipple's (1961) flux estimate is approximately 0.7, 0. 6, and 0.45, 
respectively, for a design survival probability of 0. 9 and an N g /N ratio of 0. 8, is 0. 6, 

0. 5, and 0. 4, respectively, for a design survival probability of 0. 95, and is 0. 45, 0. 4, 
and 0.3, respectively, for a survival probability of 0. 99. Though the weight of any nec- 
essary segmenting hardware is not included, substantial weight savings should be attain- 
able. 

Figure 10 shows the results when Watson's (1956) meteoroid flux estimate was used. 
For 100 segments, the minimum relative weight W /W ranged from 0. 55 to 0. 7 for a 
design survival probability of 0. 9, from 0. 48 to 0. 62 for a design survival probability of 
0. 95, and from 0. 35 to 0. 48 for a survival probability of 0. 99. While the weight savings 
of segmenting are not as pronounced as when using the more severe Whipple’s (1961) flux 
estimate, segmenting potentially can provide significant weight reduction. Again it should 
be emphasized that the manifolds were not segmented. Had the manifolds been segmented, 
even lower relative weights would have been obtained. 
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Relative radiator weight, W s /W 


I 




Ratio of surviving to total number of seg- 
ments, N s /N 


(b) Overall survival probability, 0.99. 

Figure 11. - Relative weight obtained when tube 
length was optimized, compared to original 
relative weight. Whipple's flux (1961). 


A digression at this point is in order to ex- 
plain that the segmented radiator was not opti- 
mized with respect to the tube length. Even 
though the tube length associated with the mini- 
mum radiator weight, for any given size power 
system, was not the same for each design sur- 
vival probability and meteoroid flux, the refer- 
ence radiator was considered to have tube lengths 
of 15 feet for the 1 -megawatt system, 20 feet for 
the 2-megawatt system, 30 feet for the 
5-megawatt system, and 40 feet for the 
10-megawatt system. These lengths were se- 
lected from a plot similar to figure 5 (p. 9) but 
for a design survival probability of 0. 9. The 
same tube lengths were used for all the segment- 
ing calculations. The complete optimization, 
however, where tube length was also varied was 
done for comparison by using Whipple's (1961) 
meteoroid flux for 1- and 10-megawatt-electric 
power levels and for 10 and 100 segments at each 
power level. This was done for overall survival 
probabilities of 0. 9 and 0. 99. The lightest seg- 
mented weight at each N /N ratio was used in 
conjunction with the lightest unsegmented 

weight. The new values of W /W were then 

s u 

superimposed on the original plots of figure 9 
where S was 0. 9 (N = 10, 100) and 0. 99 
(N = 10, 100), and this is shown in figure 11. 

The new relative weight was the same as for the 
nonoptimized plot for S = 0. 9 and somewhat 
lower than the nonoptimized plot for S = 0. 99. 
Thus the relative weights presented in figures 9 
and 10 may be somewhat higher than if the tube 
length had been included in the segmented radia- 
tor optimization. 

Plotting the minimum weights for segmented 
tube panels and unsegmented manifolds on a 
log -log plot against power level (fig. 12, p. 20) 
resulted in a straight line as it did for the un- 
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System power output, P e , MWe 

(a) Overall survival probability, 0.9. (b) Overall survival probability, 0.95. (c) Overall survival probability, 0.99. 

Figure 12. - Effect of meteoroid flux estimate, number of segments, and system power level on segmented radiator specific weight 

segmented case, which again indicated that the weight is proportional to the power level 
raised to some exponent. For the unsegmented case, the value of the exponent b was 
about 1. 4 and was influenced slightly by the survival probability and the choice of meteo- 
roid flux. When segmenting of the tube panels was employed, b was reduced to lower 
values. When Whipple’s (1961) meteoroid flux estimate was used, the value of b for all 
three design survival probabilities was 1. 37 for 10 segments and 1. 34 for 100 segments. 
With Watson’s (1956) meteoroid estimate and the same three design survival probabilities 
the exponent b was approximately 1. 3 for 10 segments and 1.25 for 100 segments. The 
trends are in the expected direction, because increasing the number of segments or de- 
creasing the severity of the meteoroid flux estimate decreases the amount of necessary 
protective armor. If no protective armor were necessary, the weight increase with 
power increase should be almost linear, which means that b would be unity . 

CONCLUDING REMARKS 

The salient results of the study oriented principally toward radiator weight aspects 
of large megawatt size Rankine cycle space power systems are as follows: 
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1. Radiator segmenting can offer significant weight savings especially if high survival 
probabilities are demanded and/or a high meteoroid flux estimate is used. If a radiator 
design survival probability of 0. 99 were demanded and a severe meteoroid flux such as 
Whipple's (1961) were used, 100 segments would decrease the weight to less than 0.3 of 
the weight of an unsegmented radiator . For the same design survival probability and 
meteoroid flux estimate and as few as 10 segments, the weight of the segmented radiator 
was still less than half the weight of the unsegmented radiator. With a less severe meteo- 
roid flux estimate such as Watson's (1956), the weight savings, while not as spectacular, 
are still large enough to make segmentation worthy of consideration. In this study the 
tube panels were segmented but not the manifolds. If the manifolds had been segmented, 
the results would differ considerably; undoubtedly greater weight savings would have 
been possible. 

2. The ratio of the design value of surviving segments to total number of segments 
N/N that produced the minimum weight radiator was higher than the expected value of 
0. 75. Values of 0. 8 to 0. 9 were found to yield minimum radiator weight, but this up- 
ward shift is probably due to the fact that the manifolds were unsegmented. Higher values 
of N /N imply fewer redundant segments, which results in shorter and lighter manifolds. 
Therefore, if the manifolds had been segmented, the optimum ratio of N g /N would 
probably be closer to 0.75. 

3. The radiator weight for constant cycle conditions was proportional to the gene- 

L. 

rated power raised to some exponent, W cc p“. The exponent b was a function of the 
meteoroid flux and the desired survival probability. For the unsegmented beryllium 
radiator at 1400° F, the value of b was approximately 1.4 when Whipple’s (1961) meteo- 
roid flux estimate was used, and just slightly less for Watson's (1956) meteoroid flux 
estimate. For the segmented radiator, the value of b was lower than for the unseg- 
mented radiator. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 13, 1965. 
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APPENDIX A 


SYMBOLS 


m 

o 

manifold surface area, ft 

A S 

surface area, ft 2 

A 

area vulnerable to meteoroid 

V 

2 

penetration, ft 


internal surface area of an indi- 
vidual tube, ft 2 

a s 

thin plate and spall adjustment 
in eq. (B7) 

a„ 

vulnerable area in one radiator 

V 

2 

segment in eq. (7 ), ft 

b 

exponent in eq. (1) 

C 

defined by eq. (H9) 

c 

sonic velocity in tube armor, 
ft/ sec 

D i 

internal tube diameter, ft 


manifold local internal diameter 
ft 

D o 

external tube diameter, ft 

E t 

2 

Young’s modulus, lbf/in. 

e 

element of fin 

F fin 

view factor from fin to space, 
average 

F t 

view factor from tube to space, 
average 

F > 

meteoroid flux, 
particles/ (ft 2 )(yr) 

f 

Fanning friction factor 

f(L g ) 

function of generalized length 
parameter L g 
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G 

2 

mass velocity, lbm/(ft )(sec) 

Sc 

gravitational constant, 
32.2 ft-lbm/(lbf)(sec 2 ) 

h 

fluid enthalpy, Btu/lbm 

J 

Joule's equivalent, 
778. 2 ft-lbf/Btu 

k 

thermal conductivity, 
Btu/(hr)(ft)(°F) 

L g 

generalized length parameter 
defined by eq. (C2) 

i 

distance from fin root or distance 
from manifold entrance, ft 

*fin 

half fin length, ft 

*m 

manifold length, distance along 
manifold axis between station 
of maximum diameter and 
station of minimum diameter, ft 


individual tube length, ft 

M 

total fluid mass flow through the 
radiator, lom/sec 

j( 

local fluid mass flow in individual 
manifold, lbm/sec 

m 

fluid mass flow through individual 
tube, lbm/sec 

m p 

meteoroid particle mass, g 

N 

number of radiator segments 

N P 

number of radiator panels 

% 

Reynolds number 
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N s 

design value for the number of 
surviving radiator segments 

N t 

number of tubes 

n 

constant in eq. (B7) 

P 

static pressure, psi 

P e 

system power output, kW or MW 

P(0) 

probability of receiving no 
meteoroid puncture during 
mission duration 

P 

probability that individual seg- 
ment will not be punctured by 
a meteoroid during mission 
duration 

Of in 

heat rejected per unit time from 
individual fin, Btu/hr 

Qm 

heat rejected per unit time from 
manifolds, Btu/hr 

Q T ot 

heat rejected from entire radi- 
ator, Btu/hr 

Qt 

heat extracted from working 
fluid per unit time per tube, 
Btu/hr 

%n 

heat radiated from individual fin 
per unit length per unit time, 
Btu/ (hr) (ft) 


heat radiated from tube external 
surface per unit time per tube, 
Btu/hr 

R 

local radius of manifold, ft 

S 

probability that N segments 
s 

will not be punctured by meteo- 
roids for mission duration 

s 

manifold surface arc, eq. (FI), 
ft 



temperature at midfin, °R 

T o 

fin root temperature, °R 

T w 

tube internal wall temperature, °R 

T x 

tube surface temperature, °R 

AT sc 

amount of subcooling of the con- 
densate at exit, °R 

*a 

armor thickness, ft 

£ 

difference between manifold ex- 
ternal and internal radius, ft 

^in 

fin thickness, ft 

V 

fluid velocity, ft/ sec 

•¥ 

3 

material volume, eq. (G3), ft 

V 

O 

fluid specific volume, ft /lbm 

V 

assumed average meteoroid 
velocity, ft/sec 

W 

unsegmented radiator weight, lbm 

W a 

tube armor weight, lbm 

W s 

segmented radiator weight, lbm 

W u 

unsegmented reference radiator 
weight, lbm 

X 

vapor quality 

0 

exponent of particle mass in 
eq. (4) 

y 

coefficient in eq. (B7) 

A 

increment 

€ 

surface emissivity 

V 

defined by eq. ( 6 ) 

% 

value of rj when N = N =1 

s 

e 

velocity exponent in eq. (B7) 

M 

fluid viscosity, lbm/ (ft)(sec) 
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p density, lbm/ft 8 

a Stefan-Boltzman constant, 

0. 173xl0“ 8 Btu/(hr)(ft 2 )(°R 4 ) 

r mission duration time, yr 

angles defined in appendix H 

cp density exponent in eq. (B7) 

<Pg Lockhart-Martinelli multiplier, 

eq. (D4) 

two-phase flow modulus for 
turbulent -turbulent case, 
eq. (D5) 

$ angle in eq. (H8) 

O angle in eq. (H6) 

Subscripts: 
f liquid 


fg refers to change by evaporation 

g vapor 

l midfin 

M momentum 

m manifold 

n number of times an iterative 

operation has been performed 

p meteoroid particle 

T target or protective armor 

TPF two -phase friction 

t tube 

0 initial 

Superscript: 

' refers to prime surface radia- 

tor where entire radiator is 
at one constant temperature 
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APPENDIX B 


ANALYSIS 

The calculations for this analysis were performed with a computer program, but in- 
stead of showing all the details of the program, only the basic equations that enable the 
reader to evaluate the method of analysis will be shown. Initial inputs were the radiator 
heat load, the fluid total temperature and inlet quality, the individual tube length, the 
number of panels, the desired tube pressure drop, the physical and thermal properties of 
materials, the initial number of tubes, and a guess of tube internal diameter. 

The description is for the unsegmented radiator. The differences between the equa- 
tions used for the unsegmented and the segmented radiators where the tube banks were 
segmented will be indicated at the end of this appendix. 

Before listing the equations a brief summary of the approach used may be of assist- 
ance. First, the individual tube mass flow or the total flow through the radiator divided 
by the number of tubes, the individual tube length, and the desired tube pressure drop 
were used to find the tube internal diameter. Then knowing the total tube wall area vul- 
nerable to meteoroid puncture, based on internal diameter, an armor thickness was cal- 
culated. With this first estimate for armor thickness the vulnerable area based on tube 
outer diameter was used to calculate the second estimate of armor thickness. This step 
was repeated until the ratio of two successive values was unity plus or minus a 1-percent 
tolerance. 

Next the heat radiated from the tube outer surface q^. was calculated. A surface 
temperature equal to the internal wall temperature, which was assumed to be the same 
as the fluid total temperature, was used to obtain a first estimate of the heat radiated 
from the tubes. With the first estimate of q^ the temperature drop through the armor 
was calculated. With the new surface temperature q^ was recalculated. This was re- 
peated until the surface temperature and q^ were compatible. 

Since q^. and the total heat rejection requirement were known, the heat load that had 
to be rejected by additional finning was determined. Fin dimensions were then calculated, 
and with this first estimate of fin length the manifold length and surface area were calcu- 
lated. The armor thickness was recalculated by using the total tube external surface 
area plus manifold surface area as the area vulnerable. The manifold heat rejection was 
calculated, heat radiated from the tubes recalculated, and the new fin heat load used to 
obtain new fin dimensions. This procedure was repeated until the fin geometry was such 
that the heat rejected from the tubes, manifolds, and fins to space matched the required 
radiator heat load. Then radiator weight was calculated. The total number of tubes was 
changed and the process repeated. 

The equations used in the analysis are now given. Mass flow through the individual 
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tube is 



(Bl) 


where M is total radiator mass flow and N^. the number of tubes. Since tube mass 
flow and tube length were known and a guess for tube internal diameter was made, the 
tube pressure drop was calculated by using the equations given in appendix D: 


AP t = 


AP TPF + P M 


(B2) 


which is a sum of the two -phase friction and momentum pressure changes. The AP^. 
was tested to see if 


(Desired AP^. - 0. 05 psi) < AP|. < (Desired AP^. + 0. 05 psi) (B3) 


If APj_ did not fall in the desired range, D^ was reestimated by using 





Desired AP^. 


0. 16 


(B4) 


and then AP^. was recalculated. The exponent in (B4) was an arbitrary choice based on 
some preliminary calculations. 

When the conditions of (B3) were satisfied, attention was turned to the calculation of 
the protective armor thickness required. For the first estimate it was said that the tube 
outer diameter was equal to the tube internal diameter and that manifold vulnerable area 
A m was zero: 


A v,t=V D o N t 

(B5) 

A v ~ A v, t + A m 

(B6) 


Equation (B7), from reference 11, for the armor thickness required was 
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(B7) 



where 


a s = i.75 



c 



n = 1 


v = 98 400 ft/ sec y = 2 



(The values of these parameters were taken from ref. 11 and used herein.) 
The first trial had been based on D q = D^; therefore, D q was changed to 


D 


o 


= + 2t 


a 


and equation (B7) was used to get a new t . This was repeated until 

d 



0 . 001 
12 


The heat radiated from each tube was 

\ 0 = CT€7rI Vt F t, 0 T x 

where q = 0. 75 and T x = T w were used as first estimates. Then the tube 
temperature was calculated in accordance with reference 16. 



(B8) 


(B9) 


(B10) 

surface 


(Bll) 
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When the new value for T x was used, q^ was recalculated by using equation (BIO). 
This process was repeated until 


0.999 < 



< 1 . 001 


Then the required fin heat load was calculated by using 


%in = 


( ^TOT 
*t N 



(B12) 


(B13) 


Fin thickness and fin length were then calculated by using equations (B14) and (B15) with 
a first estimate for Fjj : 




© 


f (L ) x 632. 5 2 <rk x 2e x 10 9 (— ^ 
s \1000y 


(B14) 


£ fin L g 


kt fin 


2ael0' 


9/ T 0 

aooo. 


1/2 


(B15) 


The development of equations (B14) and (B15) is shown in appendix C. With the fin length 
and D q known, the fin view factor was calculated as shown in appendix H. Then fin 
thickness and fin length were recalculated by using the new view factor. This was re- 
peated until 


(tfin) 

0.999 < — < 1.001 

( tfin )n-l 

Next the manifold surface area developed in appendix F was calculated: 


(B16) 
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A m A m, g 


A m,f = 47r£ m 


M 


a/2 


3 Wg, 


+ t_ 


+ 2^ m 


M 


a/2 


3 VrN pPfVf/ 


+ t 


(B17) 


For the unsegmented radiator the manifold vulnerable area was included with the 
tubing vulnerable area by returning to equation (B6) and recalculating everything through 
equation (B16). When the armor thickness based on the combined tube and manifold area 
had converged to within ±0. 1 percent of the previous value it was accepted. 

The heat rejected from the manifold was calculated next. The view factor was as- 


sumed to be 0. 85 and for a first estimate used T = T 

& W 

and T m,f" T w- (iT scy 2: 

Q m,g = °' 85a6A m,g T m,g 

(B18) 

^m, f = °- 85cr€A m,f T m,f 

(B19) 

O = Q + Q c 

(B20) 


where the subscript g refers to the vapor manifold and f to the liquid manifold. Then 
the manifold external temperatures were calculated by treating the manifold as a flat 
plate of area A m ^ or A m ^ and thickness t & : 


m,g 


Q t 

_ T _ mi, g a 


w 


kA. 


m ,g 


T r 

m, t 



w 


kA 


m.f 


<B21) 


(B22) 


Then Q „ and (D f were recalculated by using the calculated values for T m and 
m, g m, i J m, g 

T m f. This was repeated until 


( Q m) n 

0. 999 < — < 1. 001 

"Wn.r 

The heat to be rejected by the fin was recalculated: 


(B23) 
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(B24) 


Q TO t " Q m ‘ N t q t 


The F^. was calculated as shown in appendix H. With the new value for q^ n the pro- 
gram returned to equation (B14) to calculate a new t fin< Everything from equation (B14) 
to (B24) was repeated until 




( q fin) 

0.999 /n < 1.001 

'Wn-l 


(B25) 


Then the total vulnerable area was calculated and checked against the previous value to 
see that 


0.999 < 



(B26) 


If the condition of equation (B26) was not met, the program returned to equation (B6). 
When the ratio fell within the range specified in equation (B26), the manifold pressure 
drops and radiator weight were calculated. 

The pressure drops in the manifolds were calculated by the equations developed in 
appendix E. The fluid velocities V g and were assumed. The actual values used 
are listed in the assumptions in the text. The working fluid properties were obtained 
from reference 15: 


2 n p„vzi vn 

Ap y g, (Tg g m 

m ’S“ /«,V\1/ 


14 4g c /JiL_Y /2 

wJ 


0. 0028 + 


0.3676 


' 2M P g V g 




.0. 16 


AP. 


m,f 


144 g 


jWm 

4M V/ 2 


In. 


p irPf V f y 


0. 0028 + 


0. 3676 


^4Mp f V fl 


,0.16 


(B27) 


(B28) 
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The manifold weight equations from appendix G were 


W m,g " 2p T lrt a fi m 


4 / M \ 1/2 

3 Ivv 


+ 


W m,f _ N p fl T* t a A m 


4 / M 


a/2 


+ 


3 l N o,7rV. 


p P f ffV f y 


and the remaining weights were 


W fin " 2 p fin 1 t N t £ fin t fin 
w a = P T ‘t N t \ ( D o - < liner) 

W liner = P U „er S t N t J ( D o, liner " D i ) 
w - w m,g + w m,f + w fi„ + W a + W liner 


(B29) 

(B30) 

(B31) 

(B32) 

(B33) 

(B34) 


Modifications for Segmenting 

Segmenting was employed to reduce the amount of protective armor necessary to 
maintain the same radiator survival probability as for the unsegmented radiator. Re- 
dundant radiator sections were carried so that if an active radiator section were punctured 
by a meteoroid the damaged section would be isolated and one of the redundant sections 
activated. Thus the radiator with N segments was designed to reject the entire heat 
load from N segments, N being the design value for the number of segments which 
would remain unpunctured at the end of the mission. This requires minor modifications 
of the equations used for the unsegmented radiator. First of all, since only N g seg- 
ments out of a total of N segments are active, the number of active tubes is the frac- 
tion N /N times the total number of tubes. Thus in equations (Bl), (B13), and (B24) 
s 

where appears it is replaced by N t (N g /N). 

The other, and most important difference, is the procedure for calculating the armor 
thicknesses. The armor thickness for the manifold was calculated separately from the 
tube armor thickness. This was necessary because the manifolds were not segmented, 
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only the tube banks. The overall survival probability S that the manifolds and N out 

s 

of N independent segments would remain unpunctured is 

N 

s = P( 0 ) m P( 0 ) N = P( 0 ) V — p n (l - p) N * n (B35) 

m iN s m L^i (N - n):n: 

n=N g 

The value of S was assigned and P(0) m , the probability that the manifolds would not be 
punctured, and P(0) M , the probability that N of N segments would remain unpunc- 

JN g S 

tured, were weighted equally, or 


P(0) m = P(O) N =VS (B36) 

S 

A value of p, the probability that any individual segment would not be punctured, was 
found so that 


N 


P(0) 


N 


-£ 


N! n 


p"(l-p) N - n =V5 


s £ — * (N - n)jn. f 
n=N 


The tube vulnerable area of an individual segment is 


v N N 1 0 


(B37) 


(B38) 


Then the tube armor thickness would be calculated by using equation (B7) with this modi- 
fication. The term £o?A v y-in p(o)] 1/,3iS was replaced by the term ^afa v T/-ln pj^ 3 ^. 
Similarly the manifold armor thickness was calculated by using equation (B7) with the 
term P(0)^] 1//3 ^ replaced by [ a A m r/-ln PCO)^ 1 ^ 3 ^. Equation (B6) is de- 

leted when calculating the segmented radiator. Essentially these were the only impor- 
tant differences between the segmented and unsegmented calculation procedures. 
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APPENDIX C 


FIN LENGTH AND THICKNESS 

The equations for calculating half fin length £^ n and fin thickness, t^ n with fin 
width l held constant were taken from reference 13. One modification, adding a fin 

view factor, was made to account for the fact that the fin 
did not have an unrestricted view of space, see sketch (a). 



Substituting equation (C2) into (Cl) and simplifying yield 


Qf 


w 


in 

V F fin^t/ 


f(L )(632.46) 2 crk2el0 9 

o 


— Y 
1000 / 


(C3) 


Since has been chosen and the equivalent sink temperature was assumed to be zero, 
can be obtained from figure 8 of reference 13. With known and f(L ) obtained 
from equation (39) of reference 13, t fin can be found from equation (C3). It remains 
then to solve for £ fin in equation (C2). 
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APPENDIX D 


TWO-PHASE PRESSURE DROP IN A CONSTANT DIAMETER TUBE 

The flow regime assumed was that of a homogeneous mixture with both the liquid and 
vapor traveling at the same velocity. To obtain two-phase pressure drop, the single- 
phase multiplier of reference 17 was used. 

The tube was divided into increments and the two-phase friction static pressure drop 
for each increment was calculated by using the local vapor velocity. Because the con- 
densing process is essentially isothermal the wall heat flux was deemed to be constant. 

The energy equation as applied to the increment is 


h n " h n-l + 


V 2 - V 2 , 
n Y n-1 

2g c J 


Qt 

4— M, 

A 1 
A w 

3600 GD i 


(Dl) 


For the initial calculation the velocity term was neglected but then considered after the 
local vapor quality 


x n ~ 


h n~ h f 

h f g 


was found. The average quality for the increment is then 


x n, n-1 


x „ i + 

n-1 n 


(D2) 


(D3) 


The use of the average quality will be evident after showing the pressure drop equation. 
From reference 17 



TPF 




where cp^ 


is a function of y. 


For the assumed turbulent -turbulent case, 


(D4) 
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0.9 



(D5) 


The fluid density p and viscosity p were evaluated at inlet total temperature and pres- 
sure. The weight flow of the liquid m { and the weight flow of the vapor m g can be ex- 
pressed as a function of the vapor quality 


riij = (1 - x)m 


(D6) 


and 


in = xm 
g 


(D7) 


Thus, 


( Xtt )n,n-1 


^ 0 >f\ 0 '7 i 


0.9 


- 1 


, p f 




n, n-1 


A curve fit for the plot of <p against y+t was used to obtain cp . 

g r g 

AjL 


(AP, 


TPF 


'n, n-1 


144g c D i 


Then 


where 


(U,„-l =GX ",n-l V g 


and 


f n n-1 = 0 - 0014 + 


0. 125 




k 0. 32 


(D8) 


(D9) 


(DIO) 


(Dll) 
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Equation (Dll) is the empirical relation of Koo (ref. 16) and is valid for the range of Njj 
from 3 000 to 3 000 000. This was chosen instead of the Blasius friction factor, recom- 
mended in reference 17, because it was applicable over a wider range of Reynolds num- 
bers. Then h n was corrected for the change in velocity across the increment and the 
pressure drop for the next increment found. The friction two-phase pressure drop was 
added to the momentum pressure rise: 


A p t - -AP TpF + AP m = 



144 g c D i 


n t ° 2x 0 v fg 
144g„ 


(D12) 


n=l 
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APPENDIX E 


MANIFOLD PRESSURE DROP 


It was approximated that the fluid mass flow in the manifold decreased linearly with 
manifold length. The mass flow at any station l is then 


jt 


(1 



(El) 


where JK denotes the local fluid mass flow in the individual manifold. For the vapor 
manifold, 


J( Q 


M 

2 


and for the liquid manifold 




M 



Also 


S> i = 2R 


Secondly, it was stipulated that the manifold diameter would be varied to maintain a con- 
stant fluid velocity for the decreasing mass flow. 

The two-phase friction pressure drop in the vapor manifold is 


4£/ T pF 144g c ^. 


When Koo's empirical equation (ref. 16), is used 


f = 0. 0014 + 


0. 125 



(E2) 


(E3) 
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which is valid for 3 000 < < 3 000 000 and 


N R = 


>iVg 


Substituting equations (E3) and (E4) into (E2) yields 


'dP\ 

^di /'ppp 144g' c ^i 


2(p 2 pV 2 

Tg g 


0. 0014 +- 


0. 125 


V p 

JTg 

M 0 


0. 32 


^ 0.32 

l 


Mass flow is 




(E4) 


(E5) 


(E6) 


Rearranging (E6) and substituting (El) yield 



/ 4 J_ 

Wz, 


\l/ 2 


4 -^0 / _ i \ 

.Vg\ V 


(E7) 


Substituting equation (E7) into (E5) and assuming that the change in vapor quality is small 
enabling p , p , and cp in addition to V to be treated as constants yield equation (E5) 
as 
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Integration of equation (E8) for the vapor manifold friction pressure drop and substitution 
for yield 


AP_. 


m,g 


2 cp 2 p V 2 l 

g m 


144g 


/ 2M 

V rp g v g> 


\l/2 


0. 0028 + 


0. 3676 


For the liquid manifold, the friction drop is 


' 2M P g V g 




k 0. 16 


AP 


2p f Vjf 


m 


m,f 


144g 


4M 


a /2 


U N p PfV fi 


0. 0028 + °- 3676 


4Mp f V f , 

V “£ 2 


v0. 16 


(E9) 


(E10) 
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APPENDIX F 


MANIFOLD SURFACE AREA 


The manifold external surface area must be calculated and included as part of the 
area vulnerable to meteoroid puncture. Also, it is required to calculate the quantity of 
heat radiated to space directly from the manifold. 

Using the notation shown in sketch (b) the external surface area is 




2tt(R + t )ds 

cL 


(FI) 


ds 



di 


In equation (E7) it was shown that 



(F2) 


(F3) 
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and 

/dR\ 2 = -*0 4 _ J _\ _1 (F4) 

Vdl / 4xpVl^ \ V/ 


Substituting equations (F2) and (F4) yields equation (FI) as 




This cumbersome expression was simplified after closer examination, 
therefore f ^ » 1 and it was found that jl q/p^Y « 1; therefore, £ ^ 


First, £ m > 1, 
»J? q /4 P ttV. 


The troublesome term in equation (F6) was the inverse hyperbolic tangent. As 


^V)-° 

4„pv4/ 
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tanh" 


uT 0 


4 wpVl 


m> 


However, looking at the expression practically, if 


^° ... > 2xl0' 5 


4irpVi 


m 


then 


tanh 


-1 


[< 5.7 


1 + 


Jl 0 


4vpVl‘ 


To simplify then, the insignificant terms in equation (F6) were deleted yielding 


As- 2 ^ 


r 


fto 

3/2 

121 

2 

m 




(*n 1/2 

V npV/ m ' 


> 


(F7) 


or 


A S - 27Tl m 


i M V2 ♦ t a 

3 [irpV/ 


(F8) 


The manifold surface areas then are 


A = 4 ill 
m,g m 


M 


a/2 


+ 1„ 


3 Wg, 


(F9) 
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A m,f - 2 V*m 


2 / M 

3 ( Vp f V f/ 


a/2 


+ 1 „ 


(F10) 


A m ~ A m,g + A m, f 


(Fll) 
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APPENDIX G 


MANIFOLD WEIGHT 

As in appendix E this is the development of the general case for the individual mani- 
fold using Jt to denote local fluid mass flow in the individual manifold. The sketch in 
appendix F is also applicable here. 

The volume of material is 







dl 


(G 1) 


assuming t* = t . Simplifying and substituting equation (E7) into (Gl) yield 


r s* jr 


2t 




since the weight is 


W m = p T^ 


(G2) 


(G3) 


Substituting equation (G2) into (G3) and integrating yield 


W m = p T"Vm 


4/^0 


.1/2 


3 \fmVj 


+ t„ 


When is replaced with the proper equivalent, the manifold weights are 


(G4) 


and 


W m, g 2p T 7rt a £ m 


4 / M 


a/2 


+ 


W m,f = 


4 / M 


, 1/2 


3 N o,ttV 4 


+ t 




(G5) 


(G6) 
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APPENDIX H 


ANGLE FACTORS 


The angle factor for the heat emission from a differential element of the fin to the 
tube has been derived as a general case in reference 18. Given two surfaces, one of 
differential size and the other a cylindrical surface of arbitrary profile, whose generating 
lines are parallel to the first surface, shown in sketch (c), it was shown that the angle 
factor of the cylindrical surface A2 from the differential element dA^ is 


FdA 2 = “ ( s * n " s * n 
l 5 2 


(HI) 


For the fin and tube geometry of the radiator <£' = it/ 2 . The angle factors from the dif- 
ferential element of the fin de to the tube surfaces designated 1 and 2 in sketch (d) are 



F de, 1 


- (1 - cos <£., ) 
2 1 


and 


(H2) 


F de,2 4 <1 ' C ° S# 2> ff 13 ' 

For this study, the interradiation between the fin 
and the tube was neglected and only the radiation to 
space considered. Thus, 

F de, space = 1 ‘ F de, 1 ' F de, 2 = \ < cos *1 + cos *2> 


(H4) 
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The angle factors from 1=0 to £ = £ fin were found in increments and a mean value 
was calculated by using the trapezoidal rule: 


n=£ fin /A£ 



1 

2 


( F de, space) n l + ( F de, space 


)n _ 


A£ 


fin 


(H5) 


The angle factor from any point on the tube to the fin and adjacent tube is found similarly. 
Here the differential element (see sketch (e)) lies on the tube and the cylindrical surface 
is the fin surface plus part of the adjacent tube: 


Fj , 1 = - (1 - cos SI) 

1 2 


(H6) 


For emission directly to space, 


= 1 - F dii. i = i <* + cos n > 


: dip, space ^ dip, ± g 


(H7) 


The relation between SI and ip is shown in reference 19 to be 


cos £2 


_ C sin ;/> - 2C cos ip - (C cos ip - 1) 


(H8) 


- 2C cos ip + 1 
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where 


C 


= 21 + 



(H9) 


The mean value of the tube angle factor for 0 < ip < n/2 was obtained numerically by 
using the trapezoidal rule: 


F 


t 


n =^/2 



n=l 



77/2 



, space 



(H10) 
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